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ABSTRACT 

Context. GRB spectra appear non-thermal, but recent observations of a few bursts with Fermi GBM have confirmed previous in- 
dications from BATSE of the presence of an underlying thermal component. Photospheric emission is indeed expected when the 
relativistic outflow emerging from the central engine becomes transparent to its own radiation, with a quasi-blackbody spectrum in 
absence of additional sub-photospheric dissipation. However, its intensity strongly depends on the acceleration mechanism - thermal 
or magnetic ~ of the flow. 

Aims. We aim to compute the theiTnal and non-thermal emissions produced by an outflow with a variable Lorentz factor, where the 
power injected at the origin is partially thermal (fraction 6,1, < 1) and partially magnetic (fraction 1 - eth)- The thermal emission 
is produced at the photosphere, and the non-thermal emission in the optically thin regime. Apart from the value of 6,1,, we want to test 
how the other model parameters affect the observed ratio of the thermal to non-thermal emission. 

Methods. We followed the adiabatic cooling of the flow from the origin to the photosphere and computed the emitted radiation, which 
is a sum of modified black bodies at different temperatures. If the non-thermal emission comes from internal shocks, it is obtained 
from a multi-shell model where a fraction of the energy dissipated in shell collision is transferred to electrons and radiated via the 
synchrotron mechanism. If, conversely, the non-thermal emission originates in magnetic reconnection, the lack of any detailed theory 
for this process forced us to use a very simple parametrisation to estimate the emitted spectrum. 

Results. If the non-thermal emission is made by internal shocks, we self-consistently obtained the light curves and spectra of the 
thermal and non-thermal components for any distribution of the Lorentz factor in the flow. If the non-thermal emission results from 
magnetic reconnection we were unable to produce a light curve and could only compare the respective non-thermal and thermal spec- 
tra. In the different considered cases, we varied the model parameters to see when the thermal component in the light curve and/or 
spectrum is likely to show up or, on the contrary, to be hidden. We finally compared our results to the proposed evidence for the 
presence of a theiTnal component in GRB spectra. Focussing on GRB 090902B and GRB 10072B, we showed how these observations 
can be used to constrain the nature and acceleration mechanism of GRB outflows. 

Key words. Gamma rays bursts: general; Radiation mechanisms: thermal; Radiation mechanisms: non-thermal; Shock waves; 
Magnetic reconnection 



1. Introduction 

The first GRB spectra were obtained by the gamma and X-ray 
spectromete rs on board the IMP-6 satellite and t he Apollo 16 
spacecraft (IChne et al.lll973l: iMetzgeret al .111974 . These eariy 
observations were complemented by the large sample of 143 
spectra collected by the Konus experirn ents on the Venera probes 
from 1978 to 1980 (iMazets et al .11198 lb . It was shown that these 
spectra could be fitted by a power law with an ex ponential cut- 
off or a broken power law da me &Desail[T975l) . and various 
physical processes were invoked to expl ain this shape, such 
as optically-thin thermal bremsstrahlung (iGilman et al.lll980l) . 
Compton sca ttering of soft photons by non-thermal electrons 
dZdziarski & L amb 1986) or synchrotron emission by thermal 
or non- thermal electrons (iBrainerd & Lamb) Il987l) . Following 
the launch of the Compton Gamma-Ray Observatory, the Burst 
and Transient Source Experiment (BATSE) provided the first 
solid indications that GRBs were located at cosmological dis- 
tances and confirmed the broken power law shape of the spec- 
tra that was represented by the phenomenological Band function 
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dBand et aljl993h . Possible cyclotron or annihilation lines found 
by previous experiments were not seen by BATSE. 

In the context of the cosmological models that were devel- 
oped hereafter, the spectra were generally interpreted in terms of 
synchrotron emission from shock accelerated electrons (see e.g. 
iPiranI [199 9^ and therefore believed to be mostly non-thermal. 
The possibility that a thermal contribution cou l d also be present 
was, however, considered by iMeszaros et al ] (I2002h . and was 
supported by several observational indications. The first one 
came from the very hard low-energy spectral slopes that are 
found in some BATSE bursts during at lea st part of the evolution 
jPreece et al.lll998l;lGhirlanda et al.i2003l) . While the commonly 
observed value of the low-energy spectral index is a ~ -1, it 
reaches 0.5 to 1 in these events, suggesting the presence of a 
Rayleigh- Jeans contribution. Then Rvde (2004, 2005) proposed 
to fit all GRB spectra in the BATSE range with the combination 
of a thermal and a power law component. Using a time resolved 
analysis, he showed that during a pulse the temperature first stays 
approximately constant before decaying as a power law of tem- 
poral index close to -2/3. 

In the previous examples the identified thermal compo- 
nent represented a major contribution responsible for the peak 
of the E^N(E) spe ctrum. A different result was obtained by 
iGuiriec et al. who found a sub-dominant thermal compo- 

nent in the Ferm;-Gamma-ray Burst Monitor (GBM) spectrum 
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of GRB 100724B, accounting for a few percents of the energy 
released by the burst. The evolution of the temperature was not 
correlated to the peak energy of the non-thermal component, fit- 
ted by a Band function. Similar results have b een found in othe r 
bright GBM bursts, such as GRB 11 0721 A (iRvde et alj|2012h . 



GRB 081207, and GRB 110920 dMcGlvnn et alj l2012h . This 



improved characterisation of thermal components in the prompt 
emission of GRBs is allowed by the larger spectral coverage of 
GBM (8 keV-40 MeV) compared to BATSE (20 keV-2 MeV), 
leading to a better quality of the spectral fits. 

On the theoretical side, a thermal emission originating 
from the photosphere is a natural predictio n of most mod- 



els based on the generic fireball sce n ario (IPaczvnsl 



?sB[T986t 

iGoodmanI 119861: IShemi & PiranI Il990t iMeszaros et al.l Il993l) 
Moreover, if the acceleration of the outflow has a thermal 
origin (as would be the case if it was powered at its basis by 
neutrino-antineutrino annihilation), this photospheric emission 
would be very bright, outshining the non-thermal emission 
produc ed by internal shocks in the 1 00 keV - 1 MeV spectral 
range jDaigne & MochkovitchI 120021) . One is therefore faced 
with the following alternative: either most of the emission 
we observe is indeed this thermal component, but it has been 
Comptonised to produce a power l aw tai l in the spectrum 
at high energy (see e.g. [ Thompsonl Il994t Rees & MeszarosI 



at nign energy (see e.g. [ i nomnsoni iiyy^ iKees & Meszarosi 
l2005tlGiannios & Spruit 2007: Beloborodov 2010*) a nd comple- 
mented at low energy by additional processes (e.g. iPe'er et al] 



120061: IVurmet al.ll2011l 



origin (see e.g. 'Begelma n & Li 
20021: Iv iahakis & Konig|^^ 



or the accele ration 



2003 



has a magnetic 
Daigne & DrenkhahnI 



iKomissarov et al. 2009: 



Tchekhovskov et al. 20 lOt IKomissarov et al.ll20 10: Granot et al. 
201 1„) and the fraction eth of thermal energy in the flow is much 



smaller than unity. 

In this work we explore the consequences of this second pos- 
sibility. We self-consistently compute the photospheric thermal 
emission and the non-thermal emission from internal shocks. If 
the internal dissipation is dominated by magnetic reconnection 
rather than internal shocks, we get the non-thermal spectrum in a 
very simple parametrised way. We discuss the conditions for the 
thermal emission to show up or, on the contrary, to be hidden. 
The paper is organised as follows: in Sect. 2 we describe the ge- 
ometry and thermodynamics of the flow and explain our method 
to compute the photospheric and non-thermal emissions: our re- 
sults are presented in Sect. 3 and discussed in Sect. 4. Finally, 
Sect. 5 is the conclusion. 



2. Model description 

2.1. Geometry and thermodynamics of the flow 

We consider a schematic model where the flow that emerges 
from the central engine is accelerated by the conversion of ther- 
mal and/or magnetic energy to kinetic energy. We do not specify 
the initial geometry of the flow (which can be largely governed 
by magnetic forces), but we assume that beyond a radius 7?sph 
it becomes spherically symmetric within a cone of half opening 
angle 0. We also define the radius /?sat, where the acceleration is 
essentially complete and suppose that /?sph < ^sat- The total in- 
jected power in the flow is E, with a fraction eth in thermal form. 
The temperature Tq at the origin of the flow can then be obtained 
from 
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Fig. 1. Schematic view of the problem geometry. The flow 
emerges from the central engine through a "circular opening" 
of radius I. Beyond a radius /?sph it expands radially within a 
cone of half opening Q. The acceleration is completed at TJsat- The 
photosphere is located at R^^ and dissipation of kinetic and/or 
magnetic energy takes place at T^diss- 



where E\^o and - n( are, respectively, the isotropic injected 
power and the section of the flow at the origin (see Fig.[T}. The 
fraction of solid angle is || - ^ (we count only one jet to be 
consistent with the definition of 5o) and a is the radiation con- 
stant. We finally get 



,1/2 
-1 



pl/4 /.-1/2 



MeV, 



(2) 



with the opening angle, injected power, and radius of the jet in 
units of 0.1 rad, 10^^ erg.s"', and 10^ cm, respectively . 

We obtain the flow equations assuming that no dissipation 
takes places below the photosphere (so that the emerging spec- 
trum at transparency will be thermal only). Mass and entropy 
conservation then lead to 



pYpS = Cst, 
T 



Cst, 



(3) 
(4) 



where p is the comoving density, S the surface perpendicular to 
the flow,yS = y/c, and r = (1 -yS^-j-i/a From Eqs. [3Pl we get 



PVT^S = Cst . 



(5) 



Using Eq.|5]we can obtain the temperature at any radius R > /?sph 
even if we ignore the details of the geometry from the basis of 
the flow up to 7?sph- With S (R) = n9~ R^ and assuming that /? ~ 1 
akeady close to the origin we have 



r(/?)-rox(r2/3/;-2/3^2/3r-i/3) 

so that, in the observer frame 
TTiR) 



I +z 1 +z 
where z is the burst redshift. 



^g-2/3^-2/3^2/3p2/3^ , 



(6) 



(7) 
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At the photospheric radius R^h (supposed to lie beyond 7?sph) 
the thermal luminosity is given by 

Lth = r^, a T\R,Q c X S (7?ph) = x {0-^"R-^l"f" T^^^) . (8) 

Equations [T] an d |8] correspond to the u s ual scaling of the fire- 
ball scenario (Mdszaros & Rees| 120001: iMeszaros etaTI 120021: 
[Oaigne & Mochkovitch 2002) with, however, a modified nor- 
malization that takes into account both the geometry and the 
mixed energy content of the outflow. Convers ely, including sub- 
photospheric dissipation as in lGiannio would change the 

scaling. 

To estimate the photospheric radius we assume that most 
of the acceleration is completed at /^ptd- This is, for exam- 
ple, th e case in the simulations made by iTchekhovskov et aP 
(1201 Ol) . where the Lorentz factor sharply increases beyond the 
stellar radius, when the flow suddenly becomes unconfined. 
Then, in a first approxi mation, the photospheric radius of a 
given shell writes (e.g. iPiranI 119991: iMeszdros & ReesI l2000t 
iDaigne & Mochkovitchl2002h 



R 



ph 



kM 
8;rcr2 



2.9x10 —r cm, 

(i + cT)rl 



(9) 



where k (/cq 2 in units of 0.2 cm^.g"') is the material opacity and F 
in units of 100) the Lorentz factor of the shell. The flow keeps 
a magnetisation cr at the end of acceleration so that £■/(!+ cr) is 
the injected kinetic power ^k- In the case of a passive magnetic 
field that is carried by the out flow without contributing to its 
acceleration dSpruit et al.ll2001l) . the magnetisation cr equals 



Cpassive — (1 ~ f th) / 6th ? 



(10) 



corresponding to a pure and complete thermal acceleration. 
Efficient magnetic acceleration leads to cr < o-passivg, whereas 
o" > CTpassive corrcspouds to an inefficient magnetic acceleration, 
for instance with no conversion of magnetic into kinetic energy 
and some conversion of thermal into magnetic energy. 

When the shell reaches the photospheric radius, it releases its 
thermal energy content while a fraction of the remaining energy 
(kinetic or magnetic) can be dissipated farther away at a radius 
^diss by internal shocks for cr ^ 0.1 - 1, or reconnection for 
higher magnetisation, contributing to the non-thermal emission 
of the burst. It is therefore expected, on theoretical grounds, that 
thermal and non- thermal components both contribute to the ob- 
serve d emission (IMeszaros et al.ll2"002t iDaigne & MochkovitchI 
Hop). As mentioned in the introduction, this was already sup- 
ported by BATSE results, with new evidence now coming from 
Fermi dGuiriec et al.1 120111: IZhang et al.1 120111) . In Sect. [3] we 
present synthetic bursts showing both contributions. We explain 
below our method to compute the thermal and non-thermal emis- 
sion. 



2.2. Thermal emission 

The thermal emission can be computed from Eqs. [T] |2] El El 
and|9l for a given set of central engine parameters eth, cr, (, 6, 
Eiso, and a distribution of the Lorentz factor in the flow. Both the 
thermal luminosity and observed temperature are related to the 
injected power and temperature at the origin of the flow via the 
same factor 



<& = (r2/3/;-2/3^2/3p2m^^-2/3 pi 
V pn ph / ISO I 



«/3 
ph ■ 



(11) 



If a constant Ei^o is assumed, the luminosity and temperature 
directly trace the distribution of the Lorentz factor r(s), where 
the Lagrangian coordinate s is the distance to the front of the 
flow at the end of the acceleration stage (s/c is the ejection 
time of the shell). The expanding flow becomes progressively 
transparent (starting from the front) and the contribution of a 
shell located at s is approximately received at an observer time 
fobs = (1 + z)s/c (Daigne & Mochkovitch 2002). However, since 
the different parts of the flow do not become transparent at the 
same radius (because R^h T"-^), additional differences in ar- 
rival time of the order of 



Afobs = (1 +z) 



R 



ph 



2cr2 



^49(1 +z) 



KO.2 £'is. 



iso,53 



(l+cr)r 



(12) 



should be included. They are, however, negligible as long as 
Afobs < fvar, the typical variability time scale of the Lorentz fac- 
tor 

The value of Afobs also gives the time scale of the luminosity 
decline after the last shell of the flow (emitted by the source 
at a time r) has reached the transparency radius (high-latitude 
emission). Unless the Lorentz factor of this shell is small or the 
burst has a very short duration, the drop in luminosity for fobs > 
( 1 +z)t is ver y steep, having initia lly a temporal decay index (see 



(i -i-z)Tis ver y steep, navmg mitia ii 
e.g. Sect. 6 inl Beloborodovl201 ih 



dLog Lth 

Log fobs 



Afobs 



(13) 



See Appendix A for a short discussion of the case where Rs^ > Rph- 



This shows that in models where the prompt emission comes 
from a Comptonised photosphere, the early decay of index a ~ 
3-5 observed in X-rays cannot be explained by the high latitude 
emission and shou ld instead be related to an effective decline of 
the central engine (iHascoet et al.ll2012h . 

The expected count rate in a given spectral range and the 
resulting spectrum are obtained from the luminosity and tem- 
perature evolution (Eqs. |7] and |8). However, we do not use a 
true Planck function for the elementary spe ctrum correspond- 
ing t o a given temperat ure. As discussed in iGoodmaiil (Il986h 
and Beloborodov (2010), geometrical effects at the photosphere 
lead to a low-energy spectral index close to a = +04 instead of 
a - +1 for a Raleigh- Jeans spectrum (see also lPe'erll2008h . We 
therefore adopt a "modified Planck function" having the modi- 
fied spectral slope at low energy, an exponential cutoff at high 
energy, peaking at ^ 3.9 x A:r as a Planck function in vFy, and 
carrying the same total energy. 

2.3. Non-thermai emission 

We first estimate the non-thermal emission assuming tha t 
it comes from internal shocks (iRees & MeszarosI |1994|) . 
For a given distribution of the Lorentz factor, we ob- 
tain hght curves and spectra using the simplified model of 
IDaigne & Mochkovitch! (1998) . where the outflow is represented 
by a large nu mber of shells that inte ract by direct collisions 
only (see also 'Kobavashi e t alJ 11997 ). The elementary spec- 
trum for each collision is a broken power law with the break at 
the synchrotron energy. The adopted values for the two spec- 
tral indices at low and high energy respectively are a- - -1 
and /3 = -2.25. The expected value for a i n the fast cool- 
ing regime should normally be -1.5 (see e.g. Sari et al1ll998l: 
Ghisellini et al. 200Q), but detailed radiative models including 
inverse Compton scatt ering in Klein-Nishina regime tend to pro- 
duce harder a slopes dPerishev et al.ll200Tl: iBosniak et al.ll2009t 
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Fig. 2. Thermal and non-thermal emission from a variable outflow - internal shock framework. Top left: initial distribution 
of the Lorentz factor in the flow. Top right: thermal (red), non-thermal (blue), and total (black) photon flux in the 8 keV - 40 
MeV spectral range. Bottom left: thermal (red), non-thermal (blue), and total (black) time-integrated spectra. Bottom right: instant 
temperature (red) of the photospheric emission and instant peak energy (blue) of the internal shock emission. The dotted and dashed 
lines correspond to the temperature and peak energy averaged over time intervals of 2 and 4 s, respectively. The adopted flow 
parameters are Eiso - 10^^^ erg.s ', eth - 0.03, <t - 0.1, £ - 3 10^ cm, and - 0.1 rad; a redshift z = 1 is assumed. 



iNakaretal .1120091: iDaiene et afll201 ih. close to the typical ob- 
served value a ^ -I (Preece et al."2000'; Kaneko et a 
iNava et al.ll201 iblGolds tein et al. 2012). The global efficiency of 
internal shocks is given by the product 



/is = Q X /di 



(14) 



where /diss is the efficiency for the dissipation of energy in 
shocks and q the fraction of the dissipated energy transferred 
to electrons and eventually radiated. Typical values of /s do not 
exceed a few percents. We used this approach to compute non- 



thermal light curves and spectra in Sect. 3.1 for the case of a low 
magnetisation cr ^ 0. 1 - 1 . 

The presence of magnetic fields reduces shock efficiency and 
may e ven prev ent shock formation for cr ^ 1 (Mimic a & Alo^ 
I2OIOI: iNaravan et aLlilPlTh . Then for cr ^ 1, the magnetic field 
cannot be ignore d, and energy mu st be extract ed by magnetic 
reconnection (e.g. lThompsonlll994t Spruit et al. j2001l) . possibly 
triggered by internal shocks ("Zhang & YaniGOlTI) . Our limited 
understanding of the relevant processes does not allow a reliable 
description of the resulting emission. Considering these difficul- 
ties, we have adopted in Sect. 3.2 a very basic and simple point 



4 



R. Hascoet et al.: Prompt thermal emission in Gamma-Ray Bursts 



of view. We do not try to predict the burst profile and obtain 
the spectrum in the following way: we suppose that a fraction 
/Nth of the total injected energy eventually goes into non-thermal 
emission with a spectrum represented by a Band function with 
low- and high-energy spectral indices a - -I and P - -2.25, 
and a peak energy obtained from the Amati relation (lAmati et alJ 
I2002h 



£p ^ 130 



fmh Eis, 



1052 erg 



keV . 



(15) 



where the expon ent and normalisation values are taken from 
iNava et al] (|2012|) . The validity of the Amati relation is strongly 
debated (see e .g. iNakar & Piranl 1200 5^. 'Band & P reec3 [20051: 
lKocevskill20 1 21: ICollazzi et al.ll2012t" ]Ghirlanda et al. ll2012l) . It is 
not clear if it corresponds to an intrinsic property of gamma-ray 
bursts or if the relation results from a complex chain of selection 
effects (threshold for burst detection, and various conditions for 
the measure of the redshift and peak energy). For the purpose of 
the present study we do not address this issue and use Eq.fTSlsim- 
ply because it is approximately satisfied by the sample of long 
bursts for which the peak energy and isotropic radiated energy 
have been measured. 

3. Results 

3. 1 . Non-thermal emission from internal shocks 

To study the relative intensities of the thermal and non-thermal 
emission components we have considered as an example a flow 
ejected by a central source active for a duration itot/c - 10 s, 
where the Lorentz factor takes the form 



r(i) = 333 <{ 1 H- - cos 



SttU 

■Stot 



X exp 



s 

2sto 



(16) 



This distribution is arbitrary (see other possible examples in 
iDaigne & Mochkovitchl lT9 98|) and was adopted simply to pro- 
duce a burst made of three pulses, i.e. not too simple and not too 
complex. It is is shown in Fig.|2](top left panel) together with the 
non-thermal light curve (between 8 keV and 40 MeV; top right 
panel) resulting from internal shocks. The related photospheric 
emission is shown in the same energy range for eth = 0.03. We 
adopt a constant Eiso = 10^^ erg.s"', cr = 0.1,^ = 310* cm and 
= 0.1 rad (i.e. 1/6 = 300 km). We also represent in the bottom 
left and bottom right panels the spectrum (thermal, non-thermal, 
and global) and the temporal evolution of the instantaneous peak 
energy and temperature of the non-thermal and thermal emis- 
sions. As the thermal spectrum is a superposition of elementary 
modified Planck functions at different temperatures, its average 
spectral slope a (with A^(£') oc E") below the peak is close to -1. 
The asymptotic value a = -1-0.4 is recovered only below a few 
keV, which corresponds to the photospheric contribution with 
the lowest temperature. 

It can be seen in Fig. |2] (top right panel) that the emission is 
initially only thermal as it takes a time 



(17) 



for the first signal from internal shocks to arrive at the observer 
In Eq.[T7]the subscript "0" refers to the radius and Lorentz fac- 
tor of the first shocked shell that contributes to the non-thermal 
emission. At late times the situation is just the opposite: the pho- 
tospheric emission abruptly stops at fobs = 20 s, while the emis- 
sion from late internal shocks (both on- and off-axis) still con- 
tribute for about 10 s (observer frame). The spectrum (Fig. 2] 



bottom left) is the sum of the thermal and non-thermal contribu- 
tions. They peak at 200/(1 + z) and 1000/(1 + z) keV, respec- 
tively. Finally the plot of the temperature and of the peak en- 
ergy of the non-thermal spectrum as a function of observer time 
shows (Fig.|2] bottom right panel) that the former is more sensi- 
tive than the latter to the fluctuations of the Lorentz factor (since 

To,. - r«/3). 

We now check how these results change when we vary the 
model parameters. We divide these parameters into three groups 
describing respectively the geometry (6, /); the energetics and 
flow-acceleration mechanism (Eiso £th, cr); and the ejecta struc- 
ture {stoi, C, f ), where C and f are the contrast (C = Fmax/rmin) 
and average of the Lorentz factor distribution. 

The geometry only affects the thermal emission (for a fixed 
Eifio)- The dependence of Tobs and Lth, iso on 6 and { is weak. From 
Eqs.HEl andlwe get T^hs Q-^l^i^l'' and Lth.iso °^ r^/^^^/s 
Similarly, changing the magnetisation has a moderate effect on 
the previous results, as long as cr < 1 . As explained above, for 
cr S 1 the whole theoretical framework adopted to compute the 
non-thermal emission probably becomes invalid. 

The consequence of increasing or decreasing the thermal 
fraction eth is illustrated in Fig.[3](left panels). With eth = 0.3 the 
thermal spectrum overtakes the non-thermal spectrum between 
20/(1 + z) and 500/(1 + z) keV, and the thermal component rep- 
resents about one third of the total in the light curve. Conversely, 
with eth - 0.003 the contribution of the thermal component to 
the global spectrum and the light curve is barely visible. 

The dependence of the temperature and thermal luminosity 
on the injected isotropic power fiiso and average Lorentz factor 
r can be also obtained from Eqs.|2]l3 and|8]yielding 



Fobs cc £.-5/12 r**/\ 



-'th, iso 



(18) 
(19) 



For the non-thermal component, the relation of the peak 
energy and luminosity to the model parameters can be found 
using the simplest possible description of internal shocks that 
o nly considers the interaction of two shells of equal mass 
^Barraud et al.ll2005l) . One gets 



E,^El[^<p(c)r-^C 



■ /is ^is 



(20) 
(21) 



where fvar is of the order of one pulse duration, (fiC) depends on 
C only, and /s, defined by Eq.[T4l is given by 



/is 



e„ X 



1 + C-2VC 
1 +C 



£<■ X ediss(C) . 



(22) 



It can be seen that the thermal and non-thermal components 
behave quite diff'erently when the injected power and average 
Lorentz factor are changed. This is illustrated in Fig. |3]where we 
increase or decrease Ei^o and F by respective factors of 10 and 2, 
compared to the reference case shown in Fig. |2] 

It appears that for a given value of eth (fth = 0.03 in Fig. [3] 
middle and right panels) the thermal component becomes more 
visible when Eiso is decreased and f increased. This is a direct 
consequence of Eqs. [18] and [19] above, which can be made even 
more explicit by defining the global thermal efficiency (assum- 
ing K = 0.2 cm^ g"') 



/h = :^.4.9xl0-^eth(l+cr)2^^ 



-2/3 ,2/3 p-2/3 f.8/3 .jn^ 
-I h '^iso,53^2 '^^•'^ 
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Fig. 3. Thermal and non-thermal emission from a variable outflow - internal shock framework: impact of the outflow 
parameters. Starting from the reference case shown in Fig. |2](with a redshift z - 1) we first change the thermal fraction eth to 0.3 
(top left panel) and 0.003 (bottom left); then the isotropic power Z^iso to 10^^ erg s"' (top middle) and lO''^ erg s"' (bottom middle); 
and finally the Lorentz factor is multiplied (top right) and divided (bottom right) by a factor of 2. Each time we vary a parameter, 
all the others keep the value corresponding to the reference case. 



to be compared to the non-thermal efficiency /nih - /is approx- 
imated by Eq. [22] While changing ^iso and f affects only the 
thermal efficiency, the opposite is true for the contrast in Lorentz 
factor C. Reducing C makes internal shocks much less efficient 
and considerably softens the emitted non-thermal spectrum. In 
the limit where C — > 1 (and moreover if the Lorentz factor is in- 
creasing outwards in the ejecta) there will be no internal shocks 
and the emission will only be thermal in the absence of an alter- 
native dissipation process. 



3.2. Non-thermal emission from magnetic dissipation 

If the non-thermal emission is produced by reconnection in 
a magnetised outflow, the problem becomes very difficult, 
with no simple way to accurately follo w the process in 
time and compute a light curve (see e.g. ISpruit et al.l l2001t 



y3 



Lvutikov & BlandfordI 2003: Giannios 2008; Zhang & Yai 
201 It iMcKinnev & Uzdens kv 2012). As explained in Sect. 2.3 
above, we adopted a very simple assumption to obtain the non- 
thermal spectrum: a Band function carrying a fraction /nii, of 
the injected energy with the peak of E^N(E) obtained from the 
Amati relation. Several examples of the thermal and non-thermal 
spectra are represented in the left panel of Fig. 2] with the ther- 
mal component still being computed with the distribution of 
Lorentz factor given by Eq.[T6] The results are shown for several 
values of eth and cr, and a fixed value of the isotropic magnetic 



power at the photosphere crE/(l + cr) = lO^-'ergs"'. As ex- 
pected, the detection of the photospheric component in the spec- 
trum is favored by a high eth, a low and a high F. The spe- 
cific dependency on the magnetisation is discussed in the next 
section. 



4. Discussion 

4.1. Relative intensity of the thermal component in GRBs 

Depending on the mechanism responsible for the accelera- 
tion of the outflow in GRBs, the consequences regarding the 
photopsheric emission are very different. In a pure fireball 
(eth - 1) powered by neutrino-antineutrino annihilation (see e.g. 
Pophametal. 1999; Zalamea & Beloborodov 2011 ]), the pre- 
dicted thermal emission is very bright (Daigne & M ochkovitchI 
I2OO2I) . To agree with the current data, the spectrum should 
then be Comptonised at high energy (as a r esult of some dis 
sipative process below the photosphere, e.g. Rees & Meszaroj 
2005; Pe'eretal. 2006; Giannios 2008; Beloborodov 201 
iLazzati & Begelmanl2010l) to produce a power law tail and soft 
ened at low energy to decrease the spectral index a from a posi- 
tive value to a negative one (possibly by the presence of an addi- 
tional non-thermal component, e.g. lPe'er et al.l2006tlVunn et alJ 
201 1). An alternative is to suppose that the flow is initially mag- 
netically dominated (eth ~ 0.1). A thermal component is still 
expected to be released at the photosphere, but it will now be 
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Fig. 4. Thermal and non-thermal emission from a variable outflow - magnetic reconnection framework. In each panel we 
show a sequence of thermal (red) and non-thermal (blue), spectra (the source redshift is z = 1). The Lorentz factor distribution 
adopted for the calculation of the thermal emission is the same as in Fig. |2] (left panel) or the same but with F divided by 2 (right 
panel). The thermal emission is computed using the formalism developed in Sect. 12.21 using five values of eth^ 0.01, 0.03, 0.10, 
0.30, and 0.5. The non-thermal spectrum is simply parametrised by a Band function, with the peak energy being given by the Amati 
relation (see text) and using five values of the efficiency of the magnetic reconnection, /i-ec - /Nth(l + cr)/o": 0.01, 0.03, 0.10, 
0.30, and 0.5. In all cases, the isotropic magnetic power at the photosphere is fixed to crEiso/(^ + cr) = 10^^ erg s"'. Finally, the 
magnetisation cr at large distance is either cr = 1 (solid lines) or 10 (dashed lines). 



sub-dominant compared to non-thermal processes such as inter- 
nal shocks or magnetic reconnection. 

We have explored this second possibility in the present paper, 
making the following assumptions: (/) we supposed that the flow 
evolves adiabatically from the origin to the photosphere, i.e. we 
did not include possible sources of heating below Rp^; (//) if the 
remaining magnetisation cr at the end of acceleration is weak and 
does not prevent the formation of internal shocks, we computed 
their contribution to the emitted radiation as if cr = 0; (///) when 
cr > 1 we limited ourselves to a very simple parametrised study 
where we assumed that a fraction /nih of the injected power goes 
into the non-thermal component. 

Regarding the light curve and spectrum of the thermal emis- 
sion we obtained the following results: 

- Both the photosphere luminosity and temperature depend on 
the same factor <]? given by Eq.[TT] which directly traces the 
evolution of the Lorentz factor if the injected power stays 
constant. 

- The duration of the photospheric emission corresponds to the 
duration t of production of the relativistic wind. For t > t the 
luminosity drops rapidly on a time scale At ~ Rph/lcF^^ ~ a 
few ms for typical values of the burst parameters. 

- The global spectrum of the thermal emission is a compos- 
ite of many elementary contributions at different tempera- 
tures. Before asymptotically reaching a slope a = -1-0.4 at 
low energy it can be much softer below the peak as shown in 
Fig. 121 A time resolved spectrum will resemble more closely 
the "modified Planck function" adopted for each elementary 
collision. 



If the non-thermal emission comes from internal shocks, its 
light curve and spectrum have been computed using the simpli- 
fied approach described in Sect. 12.31 We have varied several of 
the model parameters: fraction eth of thermal energy at the ori- 
gin of the flow, isotropic injected power Eiso, and average F to 
see under which conditions the thermal component would ap- 
pear in the observed spectrum, in the internal shock scenario for 
cr < 0.1 - 1 (Fig. [3]) or in the magnetic reconnection scenario 
for cr ^ 1 (Fig. |4]i. In the internal shock framework, assuming 
K = 0.2 cm^ ■ g the ratio Q of the thermal to non-thermal 
efficiency is given by 



/is 



6.,(l+cr)^/^rf /f £,^;^3ff 
Q edi,ss(C) 



(24) 



As illustrated by this formula and in Fig. |3] increasing eth, but 
also increasing f or reducing Ei^o or C will make the thermal 
component more visible. The magnetisation cr has a low impact 
as 1 + cr ^ 1 in this scenario. In the examples shown in Fig. [3] 
the model parameters are taken from the reference case used in 
Fig. |2]and equal h = 0.3, = 1, q = 1/3, and cr = 0.1. 
The effective constrast corresponding to the initial distribution 
of the Lorentz factor plotted in the upper left panel of Fig. |2]is 
C ^ 2.5. Then Eq. [pleads to Q ^ 7.3 x 10-^eihE7^l%Tl'\ 
in reasonable agreement with Fig. |3] for instance, the upper left 
panel corresponds to eth = 0.3, £150,53 = 1, f2 - 3, and Q ^ 0.4, 
and the bottom right panel corresponds to eth = 0.03, £[^0,53 = 1, 
Fa - 1.5, and Q ^ 0.006. Note that in the E^NiE) spectrum, 
the ratio of the maxima of the two components is expected to be 
slightly higher than Q because the photospheric component has 
a narrower spectrum than the non-thermal one. 
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In the case of magnetic reconnection, the non-thermal emis- 
sion is simply parametrised by its global efficiency /nih, and we 
have 

e - 0.49 6,h ( 1 + erf I' ' /f E7^^,, rf f^l_^ , (25) 

in good agreement with Fig. |4] (/Nth,-2 being the non-thermal 
efficiency in %). Especially, Eq. |25] shows that the ratio Q in- 
creases with the magnetisation cr, which seems counter-intuitive. 
For a fixed value of E, increasing cr reduces and therefore 
decreases the photospheric radius R^^ (see Eq. |9]l. Therefore, 
for a given value of eth, the luminosity and temperature of 
the photosphere increase. However, for a given acceleration 
mechanism, one would expect an increase of cr to be associated 
with a decrease of eth, which may affect the dependency of the 
ratio Q on the magnetisation cr. For instance, in the case of a 
passive magnetic field, including cr - o-passiyg in Eq.|25] leads to 
Q - 0.49 (1 + ezf if E7^l% ff f^l_^, i.e. a decreasing 
ratio for an increasing magnetisation. Since increasing the final 
magnetisation cr tends to decrease the photospheric radius, 
it should also be noted that the acceleration of the flow may 
well be incomplete at the photosphere in high cr scenarios. As 
discussed in appendix lAl this will also reduce the photospheric 
emission. 



4.2. Comparison to observations 

We now check how these results compare to the various obser- 
vational indications of the presence of a thermal component in 
GRB spectra. 

One first indirect indication comes from the very hard 
spectral slopes a > that are so metimes observed during 
burst evolution (e .g. iGhirlanda et al.ll2003l: iBosniak et al.ll200^ 
lAbdo et aP l2009l) . Indeed, our results allow the thermal over 
non-thermal ratio (Eq. l24l i to vary with time. A locally large F 
will boost the thermal component while a low contrast C will 
reduce the non-thermal component (under the condition that the 
non-thermal emission comes from internal shocks). This can 
explain an erratic behaviour of the a slope, but observations 
often show a regular shift of a from positive to negative values 
during a single pulse. As explained in Sect. 3.1 a thermal start 
and a non-thermal ending are predicted by our models, but 
the fraction of time during which the emission is thermal is 
generally smaller than observed in the few bursts where the 
measured a slope is positive. It remains possible to adopt a 
distribution of the Lorentz factor that would extend the duration 
of the thermal emission but smoothly connecting the thermal 
and non-thermal components may not be easy. 

A second indication comes from the recent detection of pos- 
sible thermal components in a few Fermi GRBs. GRB 090902B 
is a very peculiar case of a burst showing a spectrum best fit- 
ted by a Band function with a hard low-energy photon index, 
or a multi-colour black body, together with an additional sub- 
dominant power law (Abdoetal. 2009; Ryde et al. 2010). A 
first possibility is to interpret this burst as photospheric emis- 
sion with additional dissipative processes affecting the spectrum 
by Comptonisation and possibly other non-thermal processes. 
Conversely, using the formalism presented here, we can also de- 
rive constraints on the burst parameters for alternative scenarios 
where the Band component is interpreted as the thermal emis- 
sion produced at the photosphere without sub-photospheric dis- 



sipation, and the power law as non-thermal emission produced 
above the photosphere in the optically thin regime. From the 
analysis made by Abdo et al. (2009) and Pe'er et al. (2012), the 
implied thermal luminosity and temperature are large, Lth - 
4.6 X 10^^ erg. s"' and Tobs - 168 keV. The thermal luminos- 
ity represents 53% of the total luminosity. Then, for a given eth, 
cr, and /Nth, one can deduce from Eqs.|7]|8] and|9]the isotropic 
power Z^iso, the size Rq = £/6 of the region at the base of the 
outflow (see Fig. |5] top-left panel), and the Lorentz factor F 
(see Fig.|5] botto m-left panel ) . This is a similar approach to the 
one proposed by iPe'er et ak l (l2007l) . but within the more gen- 
eral framework defined in this paper, which allows us to con- 
sider several scenarios for the acceleration of the outflow. As 
GRB 090902B is very bright, this leads to huge values of Z^iso, 
the minimum being obtained for eth - 1, which would make 
GRB 090902B a peculiar burst associated to a situation close 
to a pure fireball. For instance, for eth = i, h - 0.3, 6-i = 1 
(i.e. Rq = 300 km), and cr <s; 1 we get an isotropic power 
Eiso ~ 2.6 X lO^'* erg.s"' and a Lorentz factor F ~ 1 160. The effi- 
ciency of the photospheric emission in this case is /th - 0.18 and 
the efficiency of the non-thermal emission above the photosphere 
is /Nth - 0.16, marginally compatible with internal shocks. This 
case is labelled as F in Fig. |5](le ft column) a nd is i n good agree- 
ment with the analysis made bv lPe'er et al. 

As most GRBs do not show such a bright thermal compo- 
nent, they usually require much lower values of eth, as illustrated 
below with GRB 100724B. Then, it is worth studying the 
possibility to model GRB 090902B with eth < 1, which would 
correspond to an initially magnetised outflow. Reducing eth in 
this burst has several consequences. As illustrated in Fig. |5](left 
column), it implies a general decrease of the efficiency, except 
if the initial size Rq - (jO is very large {Rq » 3000 km), which 
is not expected for most models of the central engine of long 
GRBs, because of their short timescale (~ 1 - 10 ms) variability. 
For instance, for eth - 10"°^ (case labelled as M,isi, where we 
keep Rq - 300 km and cr <s; 1), the thermal and non-thermal 
efficiencies are reduced to /h = 0.06 and /Nth - 0.05, leading 
to an isotropic equivalent power Z^iso = 9.4 x 10^"* erg. s"'. 
It also leads to an increase of the Lorentz factor, which was 
already quite high for the pure fireball scenario (see Fig. |5] 
bottom-left panel). Therefore, eth - 0.3 -0.5 seems a reasonable 
lower limit for GRB 090902B. Note that a decrease of eth 
compared to the standard fireball leads in this case to a reduced 
non-thermal efficiency /h in good agreement with the expected 
value for internal shocks. For higher values of eth, magnetic 
reconnection, which is supposed to have a higher efficiency, 
is a better candidate, as illustrated in Fig. |5](left column). In 
addition, for a fixed value of eth, increasing the magnetisation cr 
at the end of acceleration always reduces the constraint on the 
Lorentz factor (bottom left panel). A more detailed modelling 
of GRB 090902B would be necessary to distinguish between 
these different possibilities. 

Conversely, the results of iGuiriec et al.l (1201 lb indicating 
the presence of a sub-dominant thermal component in GRB 
100724B, representing about 4% of the total flux, point towards 
low values of eth and, therefore, a magnetic acceleration. As 
shown in Fig.|5] it is diflicult to interpret GRB 100724B within 
the standard fireball scenario (thermal acceleration) that would 
imply either /Nth > 1 or Rq < 120 km (and even Rq < 40 km 
if /Nth < 0.5 is required). This is obtained assuming a redshift 
z — \, but we checked that our conclusions are unchanged for 
larger redshifts. Then a low thermal fraction eth ^ 0.01 - 0.1 is 
required in GRB 100724B, whose spectral properties are much 
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Fig. 5. Constraints on the thermal and non-thermal emission in GRB 090902B and GRB 100724B. Top: for a given thermal 
fraction e± = 10"^, 10"'^, 10"\ 10 and 1, the radius Rq - (/d at the base of the flow is plotted as a function of the non-thermal 
efficiency The corresponding thermal efficiency /th is also shown (top jc-axis). Bottom: for a given magnetisation cr = 10"', 
10"**^, 1, 10*'^, 10', 10'^, and 10^ at the end of acceleration phase, the Lorentz factor of the flow is plotted as a function of /nih (the 
unmagnetised case cr = cannot be distinguished from the case cr = 0.1). Sets of parameters representative of the different classes 
of scenarios discussed in the paper are indicated: F (eth - 1, cr = 0) (standard fireball), M,isi (logeth - -0.5, cr = 0,) and M,is2 
(log eth = -1 .5, cr = 0) (efficient magnetic acceleration: magnetisation is low above the photosphere and the dominant non-thermal 
mechanism is internal shocks), M,reci (log eth = -0.5, cr - 10), and M,rec2 (log eth - -1.5, a - 10) (magnetised flow at large 
distance, the dominant non-thermal mechanism is magnetic reconnection). The initial radius is fixed to Rq = 300 km, a typical value 
for long GRBs. The obs ervational data (thermal flux, tempera ture, ratio of the thermal over the total flux) used for the calcula tion 
(see text) are taken from lAbdo et"al] (|2009|) : IPe'er et al.l (1201 2h for GRB 090902B (left column), and from lGuiriec et al.l(l2011h for 
GRB 100724B (right column). 



more representative of the bulk GRB population than in the 
unusual case of GRB 090902B . A similar conclusion was ob- 
tained bv lZhang" &Pe'eil (120091) in the case of the very energetic 
burst GRB 0809 16C where no bright thermal component was 
detected. Assuming a passive magnetic field below the photo- 
sphere, they obtain the constraint crpassjve ^ 15 - 20, which leads 
to the more general condition eth ~ 0.05, from Eq.fTOl 

If magnetic acceleration is common in GRBs, several sce- 
narios can be discussed. In scenarios where magnetic accelera- 
tion is efficient, implying a low magnetisation at large radius and 
a dominant role of internal shocks, the thermal fraction should 
not be much larger than a few percents to avoid an unrealistic 
efficiency /Nth. For instance in the case of GRB 100724B, for 
2=1, eth = 10"'■^ h = 0.3, 6-1 = 1 (i.e. Rq = 300 km), and 
o" «: 1 (case labelled as M,is2 in Fig. |5] right column), we get 
a non-thermal efficiency /Nth - 0.06 and a thermal efficiency 



/th ^ 2 X 10"^. The isotropic kinetic power and the Lorentz 
factor in this case are Ei^o - 5.6 x 10^^ erg.s"' and F ^ 660. 
Alternative scenarios - where the flow is still magnetised at large 
radius and magnetic reconnection is the dominant mechanism to 
produce non-thermal emission - are less constrained, because 
of the uncertainties in the underlying physics. As illustrated in 
Fig. |5] for a fixed eth, increasing cr tends to reduce F, which is 
already in the typical range of a few hundreds for cr = 0. High 
values of the non-thermal efficiency /Nth ^0.1 - 0.5 (as usu- 
ally exp ected for magnetic rec onnection, see e.g. lZhang & YanI 
201 kMcKinnev & UzdenskvUIon) also require high values of 
the thermal fraction eth ^ 0.1 - 0.3, which does not seem natural 
in such scenarios of highly magnetised outflows. 

As illustrated in Sect. [3] spectra with n on-thermal and ther- 
mal components resembling those found bv lGuiriec et al.l (1201 Ih 
are easily obtained with our model, either for a photospheric -H 
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internal shocks scenario in a case of efficient magnetic acceler- 
ation, or for a photospheric + reconnection scenario if the mag- 
netisation at large distance is still large. A potential issue for 
the internal shock scenario is the moderate variation of the tem- 
perature (within a factor of 2) found in the time-resolved ana- 
lysis. To be efficient, internal shocks require large fluctuations 
of the Lorentz factors that are even amplified in the observed 
temperature (Tobs ^ T^^^, see Fig.|2|. This may suggest that the 
non-thermal emission in GRB 100724B comes from magnetic 
reconnection. This is unfortunately difficult to test in absence of 
theoretical predictions for the spectral evolution in this case. It 
should however be noted that when the temperature drops, the 
luminosity also drops so that, in practice, the temperature can 
be determined only when it is high enough. Depending on the 
time scale for the Lorentz factor fluctuations and the temporal 
resolution of the analysis, this may artificially reduce the ampli- 
tude of the measured variations of temperature. This is illustrated 
in Fig. |2] where dotted and dashed lines show the temperature 
(and peak energy of the non-thermal spectrum) averaged over 
intervals of 2 and 4 s, respectively. It remains to be tested if this 
smoothing effect can account for GRB 100724B evolution in the 
photospheric -i- internal shocks scenario. 

5. Conclusion 

We have explored in detail GRB scenarios with two episodes of 
emission: thermal emission from the photosphere without sub- 
photospheric dissipation, and non-thermal emission from inter- 
nal dissipation above the photosphere. Our results can be used to 
interpret the data and obtain constraints on the burst parameters 
or acceleration mechanism. But one faces the difficulty arising 
from the diversity of the proposed evidence for the presence of 
a thermal component in GRB spectra. In some cases this ther- 
mal component represents a major contribution to the global 
spectrum (with additional non-thermal contributions) while in 
others it is always sub-dominant, most of the emission having 
a non-thermal origin. These different situations seem to imply 
quite different magnetic over thermal energy ratios at the origin 
of the flow. However the lack of bright thermal components in 
most GRBs clearly points towards magnetic acceleration, with 
eth ^ 0.01 in most cases, and eth - 0.01 - 0.1 in less frequent 
cases such as GRB 100724B. GRB 090902B with eth ^0.3-1 
remains an exception. 

More generally, one may wonder what would be the best 
conditions for the thermal emission to show up. Apart from 
the obvious requirement that eth should be as large as possible, 
Eq. |24]may suggest looking for events with a low Ei^o and/or a 
large average Lorentz factor This, however, supposes that these 
two quantities are independent. Having F oc Zj*^ and ^ > 1 /4 
would favor both a large Eiso and F while the opposite is true 
for q < 1/4. Finally, if the non-thermal emission comes from 
internal shocks, a pure thermal spectrum can be possible even if 
the distribution of the Lorentz factor has a low contrast or if F is 
increasing outwards. 

The observation of a burst with an unambiguous photo- 
spheric signature in its spectrum would greatly help to clar- 
ify several issues in GRB physics: (i) estimating the value of 
eth would provide insight on the acceleration mechanism of the 
flow; (//) obtaining the temperature and thermal luminosity evo- 
lution would constrain the distribution of Lorentz factor and in- 
jected power; and (///) measuring the level of temperature fluc- 
tuations with a high temporal resolution would help to discrimi- 
nate between internal shocks and magnetic reconnection for the 
non-thermal emission. Isolating the photospheric component in 



the available data is, however, not an easy task: it is generally one 
among other spectral components and possibly sub-dominant, 
and does not have a simple blackbody spectrum. 

Appendix A: Incomplete acceleration at the 
photosphere 

In the case where the flow is still accelerating at the photosphere, 
the expressions for the photospheric radius, observed tempe- 
rature, and thermal luminosity will depend both on the Lorentz 
factor at the photosphere and on its value at the end of accelera- 
tion Foo. To obtain the new expressions for R^h, Tobs, and Lth we 
write that the optical depth seen by a photon produced at Rp\i and 
leaving the flow at Rq^x is equal to unity 



kM 



87rcF2/;2 



c/r = 1 . 



(A.l) 



For simplicity we suppose in this appendix that the flow is sta- 
tionary, i.e. that M is constant and that Fco is identical for all the 
shells. Then, adopting a simple parametrisation for the Lorentz 
factor 



F = 



- ph 



for Rph<R< Rs^ 



(a > 0) 



for R > R,^ 



(A.2) 



and with M = Ei^/Foo we get (neglecting terms of the order of 

(^ph/^sat)^"^') 



kEk 



Sn(2a+l)c^r^rl^Rpi, 
This finally leads to 



<ph 



Rph,c 

and 



1 /F^ 
2a +1 iFph 



(A.3) 



(A.4) 



Z^L . i^L . (2„ + 1)2/3 /^\' ^ + If' (l^f , (A.5) 

^ obs.c i^th.c \ t» / \ "sat / 

where the index c refers to the case where the acceleration is es- 
sentially complete at the photosphere. It can be seen that an in- 
complete acceleration can substantially reduce the thermal con- 
tribution. 
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